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Abstract 

Scalar  functions  defined  on  a  topological  space  ft  are 
at  the  core  of  many  applications  such  as  shape  match¬ 
ing,  visualization  and  physical  simulations.  Topological 
persistence  is  an  approach  to  characterizing  these  func¬ 
tions.  It  measures  how  long  topological  structures  in  the 
level  sets  {x  £  ff  :  f(x)  =  c}  persist  as  c  changes. 
Recently  it  was  shown  that  the  critical  values  defining  a 
topological  structure  with  relatively  large  persistence  re¬ 
main  almost  unaffected  by  small  perturbations.  This  re¬ 
sult  suggests  that  topological  persistence  is  a  good  mea¬ 
sure  for  matching  and  comparing  scalar  functions.  We 
extend  these  results  to  critical  points  in  the  domain  by 
redefining  persistence  and  critical  points  and  replacing 
sub-level  sets  {x  €  fl  :  f(x )  <  c}  with  interval  sets 
{x  £  0  :  a  <  f(x)  <  6}.  With  these  modifications  we 
establish  a  stability  result  for  domain  points  that  can  be 
used  for  matching  two  scalar  functions. 

1  Introduction 

A  scalar  field  is  a  scalar  function  /  :  fi  — >  R  defined 
on  some  topological  space  fi.  Examples  of  scalar  fields 
are  fluid  pressure  in  computational  fluid  dynamics  sim¬ 
ulations,  temperature  in  oceanographic  or  atmospheric 
studies,  and  density  in  medical  CT  or  NMR  scans.  A 
level  set  of  a  scalar  field  is  a  set  of  points  with  the 
same  scalar  value,  i.e.,  {x  €  fl  :  f(x)  =  c}.  One 
way  of  deriving  quantitative  information  about  scalar 
fields  is  by  studying  the  topological  structures  of  its 
level  sets  or  the  regions  bounded  by  level  sets,  such  as 
{x  €  Q  :  f(x)  <  c}.  The  mathematical  field  of  Morse 
Theory  is  the  study  of  these  topological  structures. 

Among  the  most  basic  problems  on  scalar  fields  is 
simplifying  a  scalar  field  for  compact  representation, 
identifying  important  features  in  a  scalar  field,  and  char¬ 
acterizing  the  essential  structure  of  a  scalar  field.  Ex¬ 
tracting  and  representing  the  topological  structure  of  the 
level  sets  is  one  way  of  approaching  all  these  problems. 
However,  this  topological  structure  may  contain  “small” 
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topological  features  which  are  insignificant  or  caused  by 
noise.  Small  topological  features  should  be  removed 
in  simplification  and  ignored  in  characterizing  essential 
structure  or  identifying  important  features.  How  does 
one  determine  which  topological  features  are  small? 

Edelsbrunner,  Letscher  and  Zomorodian  in  [5]  intro¬ 
duced  the  notion  of  topological  persistence.  As  c  G  R 
increases,  topological  features  appear  and  disappear  in 
the  set  {x  €  fl  :  f(x)  <  c }.  If  a  topological  fea¬ 
ture  appears  at  “time”  a  and  disappears  at  “time”  b,  then 
its  persistence  is  the  difference,  b  —  a,  between  these 
two  times.  Edelsbrunner  et  al.  in  [5]  use  homology 
groups  over  Z/2Z  to  define  topological  features.  Carls- 
son  and  Zomorodian  [9]  showed  how  topological  persis¬ 
tence  could  be  computed  for  homology  groups  over  any 
fields. 

At  the  core  of  various  application  areas  such  as  shape 
matching  and  visualization  is  the  problem  of  character¬ 
izing  and  comparing  scalar  fields.  Topological  persis¬ 
tence  gives  one  approach  to  comparing  such  fields.  Two 
fields  are  similar  if  they  have  matching  topological  fea¬ 
tures  with  approximately  the  same  persistence.  This 
approach  to  comparing  fields  makes  sense  only  if  per¬ 
sistence  remains  stable  under  relatively  small  perturba¬ 
tions  of  the  scalar  fields.  Cohen-Steiner,  Edelsbrunner 
and  Harer  [3]  proved  that  “large”  persistence  values  re¬ 
main  almost  unaffected.  More  precisely,  let  scalar  field 
/  :  fl  — >  R  be  a  small  perturbation  of  field  /  :  fl  -A  R, 
(i.e.,  \f(x)—f(x)  |  <  Jforallx  6  H.)  If  /  has  a  topolog¬ 
ical  structure  with  relatively  large  persistence  which  ap¬ 
pears  at  a  and  disappears  at  b,  then  /  has  a  corresponding 
topological  structure  which  appears  around  a  and  disap¬ 
pears  around  b. 

Critical  values  are  the  range  scalar  values  where  the 
topological  structure  of  the  level  sets  changes.  Cohen- 
Steiner  et.  al.  showed  that  the  critical  values  for  struc¬ 
tures  with  large  persistence  remain  stable  under  small 
perturbations  of  the  scalar  field.  Scalar  fields  also  have 
critical  points,  points  in  the  domain  which  change  the 
topological  structure  of  the  level  sets.  It  is  natural  to  ask 
if  critical  points  for  structures  with  large  persistence  re¬ 
main  stable  under  perturbations  of  the  field.  If  two  scalar 
fields  are  close,  then  are  their  significant  critical  points 
“close”? 

In  this  paper  we  revisit  topological  persistence  and  es¬ 
tablish  a  stability  result  in  terms  of  the  critical  points  in 
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the  domain.  There  are  two  obstacles  that  we  must  over¬ 
come.  First,  if  we  look  only  at  the  topological  structure 
of  {x  €  Cl  :  f(x)  <  c},  then  there  is  no  such  stability 
result. 

Consider  Figure  1 .  Let  /  be  a  function  defined  on  the 
surface  in  K3  shown  in  the  figure  where  f(x)  is  the  21- 
coordinate  of  the  point  x.  The  set  {x  €  Cl  :  f(x)  < 
/(r)}  is  homeomorphic  to  a  cylinder  with  circles  C\  and 
C2  bounding  each  end.  The  first  homology  group  (II \)  is 
generated  by  the  circle  c\ .  (It  is  also  generated  by  circle 
C2.)  Assume  that  the  two  maxima  p  and  q  have  2- co¬ 
ordinates  that  are  almost  equal  except  that  f(p)  <  f(q). 
Since  f(p)  <  f(q),  the  first  homology  group  becomes 
trivial  in  the  set  {x  €  fl  :  f(x)  <  /(p)}.  Loosely  speak¬ 
ing,  the  cycle  generated  by  <7  at  r  gets  destroyed  at  p. 
According  to  the  framework  of  [5],  we  get  a  persistent 
value  pair  [f(r),f(p)\. 

Now  consider  a  slightly  perturbed  /  denoted  as  /.  Set 
/  equal  to  /  everywhere  except  in  the  vicinity  of  p  and 
q  where  /  is  perturbed  so  that  f(q)  <  f(p).  Let  p'  and 
q'  be  new  maxima  close  to  p  and  q  respectively  for  /. 
The  cycle  generated  by  C\  at  r  gets  destroyed  at  q1  in 
{x  £  fi  :  f(x)  <  f(q')}-  We  get  a  persistent  value  pair 
[/(r),  f(q')]  for  /.  Since  f(p)  is  close  to  /(</'),  the  two 
persistent  value  pairs,  [/(r),  f(p )]  and  [/(r),  f(q')\,  are 
close,  confirming  the  Cohen-Steiner  et  al.  result.  How¬ 
ever,  the  points  p  and  q'  are  not  close  in  any  sense  in  the 
domain. 


Figure  1:  Set  {x  :  f(x)  <  f{r)}  is  homeomorphic  to 
a  cylinder  with  ends  at  circles  C\  and  c-i.  Circle  Ci  gen¬ 
erates  the  first  homology  group  of  {x  :  f(x)  <  /(r)}. 
Since  f(p)  <  f(q ),  this  homology  group  is  destroyed 
at  pin  {x  :  f(x)  <  /(p)}-  If  /  is  a  slight  perturbation 
of  /  where  f(q')  <  f{p),  then  this  homology  group  is 
destroyed  at  q'  in  {x  :  f(x)  <  /(</')}.  Points  p  and  q1 
are  far  apart  in  the  domain. 

Instead  of  considering  only  sets  {x  £  Cl  :  f(x)  <  c} 
which  are  bounded  from  above  by  a  single  level  set, 
we  will  consider  interval  sets  that  are  bounded  from 
above  and  below  by  level  sets.  The  sets  we  use  are 
{x  £  Cl  :  a  <  f(x)  <  b}.  This  is  one  crucial  devia¬ 
tion  we  make  from  the  set  up  in  earlier  works  [3,  5,  9]. 


It  also  leads  to  slightly  different  definitions  for  critical 
points  and  persistence. 

Returning  to  Figure  1 ,  the  first  homology  group  of  the 
set  {x  :  f(r)  <  f(x)  <  /(p)}  (open  at  the  top)  has  two 
distinct  generators,  one  given  by  circle  Ci  and  one  by  cir¬ 
cle  C2-  Point  p  destroys  the  homology  group  generated 
by  ci  while  point  q  destroys  the  one  generated  by  C2- 
In  the  perturbed  held  /  where  f(q')  <  /(p'),  a  similar 
thing  happens.  Thus  p  and  p'  (also  q  and  q')  are  critical 
points  for  /  and  /  respectively  and  destroy  homology 
groups  with  approximately  the  same  persistence. 

The  second  problem  in  stability  of  critical  points  is 
that  corresponding  critical  points  can  be  arbitrarily  far 
apart.  Consider  functions  /  and  /  in  Figure  2  where 
|  f(x)  —  f(x)  |  <  6  for  all  x  £  Cl.  The  maxima,  p  and 
q,  of  /  and  /,  respectively,  can  be  made  arbitrarily  far 
apart  even  as  S  is  made  arbitrarily  small. 
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Figure  2:  The  maximum  p  for  a  real  valued  function  / 
has  moved  by  large  distance  even  for  an  arbitrarily  close 
approximant  /. 

Instead  of  using  a  metric  in  the  domain,  we  use  the 
range  to  determine  neighborhoods  of  points.  A  neigh¬ 
borhood  of  point  p  is  the  connected  component  of  {x  € 
Cl  :  f(p)  —  71  <  x  <  f(p)  +  72}  containing  p.  A  point 
which  is  in  this  neighborhood  for  small  values  of  71  and 
72  is  “close”  to  p.  Note  that  points  p  and  p'  in  Figure  2 
are  “close”  in  this  sense,  but  points  p  and  q’  in  Figure  1 
remain  far  apart.  We  show  that  if  p  destroys  a  persis¬ 
tent  homology  group  in  /,  then  the  neighborhood  of  p, 
{x  £  Cl  :  f(p)  —  71  <  x  <  f(p)  +  72},  contains  a  point 
q  which  destroys  a  persistent  homology  group  in  /.  The 
values  of  71  and  72  depend  upon  the  persistence  of  the 
homology  group  and  the  difference  <5  between  /  and  /. 

Theorem  1,  one  of  our  main  results,  states  that  ev¬ 
ery  destroying  critical  point  of  /  contains  in  its  “neigh¬ 
borhood”  a  similar  destroying  critical  point  of  /.  How¬ 
ever,  to  construct  a  matching  of  critical  points,  we  need 
each  destroying  critical  point  of  /  to  be  in  only  one  such 
neighborhood.  We  establish  this  stronger  result  for  criti¬ 
cal  points  which  are  local  maxima  of  functions  on  man¬ 
ifolds.  More  specifically,  we  show  that  the  neighbor¬ 
hoods  of  local  maxima  with  large  persistence  are  pair¬ 
wise  disjoint. 
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2  Definitions  and  assumptions 

2.1  Homology  groups 

For  a  topological  space  X,  the  fcth  homology  group 
Hk(X)  is  an  algebraic  encoding  of  the  connectivity  of 
X  in  the  Art h  dimension.  For  a  good  exposition  on  ho¬ 
mology  groups  we  refer  to  Hatcher  [6],  We  will  use  the 
singular  homology.  Although  homology  groups  are  de¬ 
fined  for  coefficients  drawn  from  any  ring,  we  will  con¬ 
sider  only  fields  such  as  M.  Q,  Zp  =  Z/pZ  for  a  prime  p 
as  in  the  previous  works  [5,  9].  As  discussed  in  Carlsson 
and  Zomorodian  [9],  computing  the  persistent  homology 
groups  over  non-fields  is  an  unsolved  and  perhaps  in¬ 
tractable  problem.  Over  fields  the  homology  groups  are 
vector  spaces  and  the  rank  of  Hk (W),  denoted  fdk(X), 
is  called  the  fcth  Betti  number  of  X. 

A  continuous  map  /:  X  — >  Y  between  two  topo¬ 
logical  spaces  X  and  Y  induces  a  homomorphism,  say 
fk,  between  their  homology  groups,  Hk (X)  4  Hk(Y). 
This  property  is  carried  over  the  composition  of  maps, 
that  is,  (/  o  g)k  =  fk  o  gk.  In  our  case,  the  maps  be¬ 
tween  spaces  will  be  inclusions  maps.  This  means,  if 
X  C  Y,  we  will  consider  the  map  Hk{X)  4  Hk(Y) 
where  l  is  induced  by  the  inclusion  map  i :  X  — >  Y . 
From  now  on,  we  take  the  liberty  of  dropping  the  sub¬ 
script  k  from  Hk(X)  when  it  is  clear  from  the  context. 
For  X  C  Y  the  relative  homology  group  of  Y  with  re¬ 
spect  to  X  is  given  by  H(Y,X )  =  H(Y)/l(H(X)) 
where  H(X)  — >  H(Y)  is  the  homomorphism  induced 
by  inclusion  t :  X  — >  Y . 

A  sequence  of  groups  CJj  connected  by  homomor- 
phisms  form  an  exact  sequence  if  any  two  consecutive 
homomorphisms  in  the  sequence 

■  ■  ■  — f  Gi  — ^  Gi-\.  i  — ^  Gi-\. 2  — >  . . . 
satisfy  the  property  that 

Im  ti  =  Ker  li+i. 

We  will  use  a  specific  type  of  sequence  called  Mayer- 
Vietoris  sequence  which  is  known  to  be  exact. 

Let  A,  B  C  X  so  that  X  is  the  union  of  the  interiors 
of  A  and  B  and  D  =  A  n  B.  The  sequence 

Hk(D)  4  Hk(A)  ®  Hk(B)  4  Hk(X)  4  Hk^(D) 

is  exact  and  is  called  the  Mayer-Vietoris  sequence  [6,  p. 
149].  The  map  d  is  the  connecting  homomorphism  given 
by  boundary  maps  [6,  p.  1 16]. 

2.2  Notation 

We  use  the  following  notation  to  define  the  region 
bounded  by  /_1(o)  and  /_1(6).  For  a,  b  £  R  and  func¬ 


tions  /  and  g.  let 

Fb  =  {x  £  12  :  a  <  f(x)  <  6}  and 

Gba  =  {x  £  0  :  a  <  g(x)  <  6}. 

In  our  results  and  proofs  we  need  the  space  Fb  and 
Gba  closed  at  the  bottom.  So,  we  define 

Fba  =  {x  £  0  :  a  <  f(x)  <  6}  and 

G ^  =  {x  £  12  :  a  <  g(x)  <  b}. 

Notice  that  a  could  be  — oo  and  b  could  be  oo. 

2.3  Destruction 

Let  ft  be  a  topological  space.  For  X  C  ft  and  Y  C  11, 
set  Y  destroys  non-zero  h  €  H(X)  if  the  image  of  h 
under  the  mapping  H(X)  — >  II (X  U  Y)  is  zero.  In 
particular,  if  q  is  a  point  in  12,  point  q  destroys  non¬ 
zero  h  £  H(X)  if  the  image  of  h  under  the  mapping 
H(X)  ->  H(X  U  {<?})  is  zero. 

If  X  C  Z  C  12  and  Y  C  12,  then  we  say  that  Y 
destroys  the  image  of  hx  £  H(X)  in  H (Z),  if  hz  £ 
H{Z)  is  the  image  of  hx  under  the  mapping  H(X)  — > 
H(Z)  and  Y  destroys  hz.  We  encounter  this  situation 
repeatedly  where  X  is  some  level  set  /  L  (a)  and  Y  is 
a  point.  For  brevity,  we  say  that  point  q  destroys  h  £ 
/_1(a)  if  point  q  destroys  the  image  of  h  in 

A  function  /  :  12  — >  M  is  point  destructible  if  when¬ 
ever  h  £  H(F_ba )  is  destroyed  by  /_1(6),  then  h  is  de¬ 
stroyed  by  some  point  q  £  /_1(6). 

2.4  Maps  and  spaces 

We  will  be  dealing  with  continuous  functions  on  a  com¬ 
pact,  connected  topological  space,  12.  We  need  some 
conditions  that  these  functions  will  be  well-behaved,  i.e. 
have  properties  similar  to  Morse  functions.  However,  we 
do  not  want  to  restrict  ourselves  to  differentiable  func¬ 
tions  or  to  Morse  functions. 

For  a  function  /  :  12  — >  R  and  a  £  R,  let 
=  {x  £  ft  :  a  —  e  <  f{x)  <  a  +  e}  denote 
the  open  e-neighborhood  of  /_1(a).  The  first  property 
we  require  is  that  the  topology  of  /_1(a)  is  similar  (iso¬ 
topic)  to  the  topology  of  a  e-neighborhood  of  /_1  (a)  for 
suitably  small  e.  The  second  property  is  that  /  is  point 
destructible.  These  properties  are  similar  to  the  Morse 
condition  that  critical  points  are  isolated.  We  define  the 
first  property  more  formally  below. 

Represent  the  unit  interval  [0, 1]  by  I.  Subspace 
X  C  F  is  a  strong  deformation  retract  of  Y  if  there 
is  a  continous  cf>  :  X  x  I  — »  X  such  that  f(y,  0)  =  y 
and  <f>(y,  1)  £  X  for  all  y  £  Y  and  <j>(x,t)  =  X  for 
all  x  £  X  and  t  £  R.  In  other  words,  <j>  continuously 
deforms  Y  into  X  without  moving  any  points  in  X.  If 
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X  is  a  strong  deformation  retract  of  Y,  then  II  (X )  is 
isomorphic  to  H (Y). 

We  say  that  the  continuous  function  /  :  LI  — »  M  is  LR 
(locally  retractible)  if  for  all  a  E  M,  there  exists  some 
eo  >  0  such  that  for  all  e  <  eo,  the  set  f^1  (a)  is  a  strong 
deformation  retract  of  Nc(f~1(a))  =  {x  :  a  —  e  < 
f(x)  <  a  +  e}. 

Piecewise  linear  functions  on  finite  simplicial  com¬ 
plexes  are  LR.  If  a  continous  function  is  Morse,  then  the 
function  is  LR.  (See  Milnor  [7,  pp.  12-20]  for  a  proof. 
Milnor  actually  proves  that  { x  :  f(x )  <  a}  is  a  de¬ 
formation  retract  of  {x  :  f(x)  <  a  +  e}  but  his  proof 
also  shows  that  /_1(a)  is  a  strong  deformation  retract 
of  {x  :  a  —  e  <  f(x)  <  a  +  e}.)  Some  properties  of 
LR  functions  will  be  used  in  our  proofs.  We  state  them 
in  Lemma  1  and  Lemma  2.  The  proofs  are  given  in  the 
appendix. 

Lemma  1  says  that  if  /  is  LR,  then  we  can  replace  {x  : 
a  <  J(x)  <  &}  by  suitably  chosen  small  neighborhoods 
without  changing  its  homology. 

Lemma  1.  If  continuous  function  f  :  Ll  — >  M  is  LR,  then 
for  every  a,  b  E  R  where  a  <  b,  there  exists  an  such 
that  for  all  c  <  eg,  set  {x  :  a  <  f{x)  <  b}  is  a  strong 
deformation  retract  of  {x  :  a  <  f(x)  <  b  +  e}  and  set 
{x  :  a  <  f{x)  <  b }  is  a  strong  deformation  retract  of 
{x  :  a  —  e  <  f(x)  <  b}. 

Let  non-zero  h  E  H(f~l  (a))  be  destroyed  by  Fff.  If 
/  is  LR,  then  h  is  destroyed  by  {x  :  a  <  f(x)  <  6}  for 
some  b  >  a.  Equivalently,  the  image  of  h  in  H(F_ba )  is 
destroyed  by  /_1(6). 

Lemma  2.  Let  f  :  f I  — >  R  he  a  continous,  LR  function. 
For  any  non-zero  h  E  H(f^1(a)),  ifH  (Fff )  destroys  h, 
then  for  some  b  >  a,  the  image  of  h  under  the  mapping 
H (/_1  (a))  — >  H (F_ba)  is  destroyed  by  f  (b). 

Lemma  2  has  the  following  corollary. 

Corollary  3.  Let  f  :  fl  — >  R  be  a  continuous, 
point  destructible,  LR  function.  For  any  non-zero  h  E 
H(f^1(a)),  if  H(Fff)  destroys  h,  then  some  q  E  Fl 
destroys  h. 

Proof.  By  Lemma  2,  there  exists  some  b  >  a  such  that 
the  image  h'  of  h  under  the  mapping  H(f~1(a))  — 
H (. F_ba)  is  destroyed  by  /_1  (6).  Since  /  is  point  destruc¬ 
tible,  there  is  some  point  q  E  /_1(6)  which  destroys  h' 
and  h.  □ 

3  Persistence 

Intuitively,  the  persistence  of  a  point  p  E  0  is  the  “age” 
of  the  “oldest”  homology  element  destroyed  by  p.  For¬ 


mally,  the  persistence  of  point  p  E  ft  is 

II  f(p)  =  limsup{/(p)  —  a:  p  destroys  some 
non-zero  h  E  H(f^1  (a))}. 

We  use  lim  sup  in  place  of  max  because  it  is  possible 
that  p  destroys  non-zero  elements  of  H (/_1  (a  +  e))  for 
any  e  >  0  but  not  elements  of  H (/_1  (a)). 

If  point  p  has  persistence  7,  and  f(p)  —  7  <  a  < 
f(p),  does  p  destroy  some  element  of  i7(/_1(a))?  As 
we  show  below,  the  answer  is  yes. 

We  will  need  to  use  subsets  of  fi  which  are  bounded 
by  two  different  functions  in  our  discussion  of  persis¬ 
tence.  The  set  X  =  {x  :  a  <  f(x)  and  g{x)  <  b} 
is  properly  bounded  by  /  and  g  at  a  and  b  respectively, 
if  the  level  sets  /_1(a)  and  5_1(6)  are  disjoint  and  X  is 
non-empty.  Notice  that  X  includes  /  1  (a)  at  the  bottom 
but  not  g-1  (b)  at  the  top. 

In  Lemma  4  we  will  establish  a  result  relating  three 
functions  and  the  spaces  delimited  by  their  level  sets. 
Later  we  will  set  these  level  sets  only  to  those  of  two 
functions  /  and  g  that  are  used  to  establish  the  stability 
result. 

Let  /1,  /2,  fs  :  Q  — >  R  be  three  continuous  LR  func¬ 
tions.  Let  q  be  a  point  in  ft  and  let  03  =  f$(q ).  Let 
X  =  {x  E  11  :  ai  <  fi (x)  and  fs{x)  <  <23}  and 
X'  =  {x  E  fi  :  a2  <  h{x)  and  /3(x)  <  be 
properly  bounded  subsets  of  fl  where  X'  C  X  and 
El  /^1(« 2)  =  0-  see  Figure  3.  The  essence 
of  the  following  lemma  is  that  if  q  destroys  a  generator 
in  H (ff1{a,i)),  then  it  necessarily  destroys  some  gener¬ 
ator  in  H (ff1  (02))-  Some  subtle  aspect  of  this  destruc¬ 
tion  is  narrated  in  the  caption  of  Figure  3. 


Figure  3:  Sets  X  =  {x  E  LI  :  a±  <  fi(x)  and 
fz(x)  <  <23}  and  X1  =  {x  E  LI  :  a 2  <  fo{x)  and 
fz(x)  <  a^}.  Cycle  C  generates  hx  E  H (X)  and 
h'x  E  H(X')  where  both  hx  and  h'x  are  destroyed  at 
q.  Elements  hx  and  hx  are  the  images  of  some  hi  E 
H{ff  1(«i))  and  some  E  H(ff1{a2)),  respectively. 
The  mapping  H(X')  -»  H (X )  sends  h'x  to  hx.  Cycle 
C'  generates  hx  E  H(X)  and  hx  E  H(X')  where  hx 
is  destroyed  at  q  but  is  not,  even  though  the  mapping 
H(X')  -»  H(X)  sends  ft"  to  hx. 
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Lemma  4.  If  q  destroys  some  non-zero  hx  G  H  (X) 
and  hx  is  the  image  of  some  h\  G  H(f~l(ai)),  then 
hx  is  the  image  of  some  non-zero  hx  G  H(X')  which 
is  destroyed  by  q.  Moreover,  h’x  is  the  image  of  some 
non-zero  h2  G  H(f~1(a2)). 

Proof.  Let  hx  €  H(X)  be  the  image  of  some  h\  € 
H(f^1(ai))  under  the  mapping  iJ(/_1(ai))  -3  H(X) 
where  hx  is  destroyed  by  q.  By  Mayer- Vietoris,  the  se¬ 
quence 

H(X’)  -3  H(X'  U  {<?})  ®  H(X)  -3  H(X  U  {<?}) 

is  exact.  Since  the  induced  mapping  H{X)  — >  H (X  U 
{</})  sends  hx  to  0,  the  mapping  H( X'  U  {g})  ® 
H(X)  — >  H(X  U  {</})  sends  (0  ®  hx)  to  0.  Since  the 
sequence  is  exact,  there  is  some  hx  whose  image  is  (0  ® 
hx)  under  the  mapping  H(X')  -3  H(X'U{q})(BH(X). 
Thus  hx  is  the  image  of  h'x  and  h'x  is  destroyed  by  q. 

We  now  prove  that  h'x  is  the  image  of  some  /i2  G 
H(ff1(a2)).  Since  /2  is  LR,  there  exists  some  ex  >  0 
such  that  H(ff1(a2))  -3  H(Ncj  (ff1(a2)))  is  an  iso¬ 
morphism  for  all  e'  <  By  Lemma  1,  H(X')  -3 
H(X'  U  Nei  (/2^1(«2)))  is  an  isomorphism  for  all  e'  < 
e2.  Since  ff1(ai)  fl  f21(a 2)  =  0-  there  is  some  €3  such 
that  Nes  (ff1  (a2))  Q  H(X).  Let  e  be  the  smaller  of  ei, 
e2  and  e3. 

Let  Y  equal  Ne(ff1(a2)).  Let  Z  =  {1  G  SI  :  «i  < 
fi  (, x )  and  /2  (x)  <  a2  +  e}  and  Z'  =  {x  G  fl  :  a2  —  e  < 
f2(x)  and  /3(a:)  <  a3}.  The  following  commutative 
diagram  gives  the  relevant  mappings  between  homology 
groups: 

h  1  G 

tf(/rV)) 

h-i  G  H(ffl(a2))  = =H{Y ) - »■  H(Z)  3  hz 


t  Mg  t 

K  G  iT(X')  JT(Z') - ^  9  hx 


must  be  in  the  image  of  some  /iy  G  H(Y).  Since 
H (ff1(a2))  -3  H (Y)  is  an  isomorphism,  there  is  some 
h2  G  H (ff1  (o2))  whose  image  is  hy  under  the  map¬ 
ping  H(Jf  \a2 ))  ->  H(Y). 

All  the  mappings  of  homology  groups  are  induced  by 
the  inclusion  mapping  and  thus  the  diagram  is  commu¬ 
tative.  Since  H  (X1)  -3  H  ( Z ')  is  an  isomorphism,  the 
inverse  mapping  H(Z')  -3  H(X')  takes  h'z  to  h'x.  The 
mapping  H{f2\a2))  -3  H{Y)  -3  H(Z')  -3  H(X') 
sends  element  h2  G  H(f2  1(a2)  to  h'x  G  H{ X').  Thus 
h'x  G  H(X')  is  the  image  of  some  non-zero  h2  G 
H{f2\a2).  □ 

Setting  fi  =  f2  =  /3  gives  the  following  corollary: 

Corollary  5.  Let  f  :  LI  — }  R  be  a  continuous,  LR  func¬ 
tion.  If  q  G  fl  has  persistence  7,  then  for  every  a  where 
f(a)  —  7  <  a  <  f(q),  point  q  destroys  some  element  of 

"'■I  '(«)). 

4  Stability 

In  this  section  we  prove  one  of  our  main  results.  Theo¬ 
rem  1.  Let  /  and  g  be  two  functions  defined  on  fl.  We 
say  |/  —  g\  <  6  if  |  f(x)  —  g(x) \  <  5  for  all  x  G  0. 
We  show  that  if  q  is  a  destructor  for  /  with  persistence 
7  >  28,  there  is  a  point  q'  which  is  a  destructor  for  g 
where  q  and  q'  lie  in  the  same  connected  component  of 
the  space  (Wherever  we  use  the  term  con¬ 

nected  component,  we  always  mean  path  connected.) 
Moreover,  the  values  f(q)  and  g(q')  are  close.  This  the¬ 
orem  not  only  relates  q  and  q'  in  the  range  as  in  Cohen- 
Steiner  et  al.  [3]  but  also  in  the  domain. 

We  will  need  the  concept  of  chains  and  cycles  that 
define  the  homology  groups.  Chains  are  formal  sums  of 
maps  from  standard  simplices  into  the  domain  fl  and  cy¬ 
cles  are  chains  which  have  no  boundary.  The  boundary 
of  a  chain  is  always  a  cycle.  For  details,  see  Hatcher  [6]. 

Lemma  6.  If  non-zero  h  G  H(F_ba)  is  destroyed  by  point 
q  and  cycle  C  C  F_ba  generates  h,  then  C  and  q  lie  in  the 
same  connected  component  of  cl(F_ba). 


Element  hx  is  the  image  of  hi  G  H(ffl(a  1)).  Let  hz  Proof.  Since  q  destroys  h,  cycle  C  is  the  boundary  of 

be  the  image  of  hi  under  the  mapping  H(f^1(ai))  -3  some  chain  D  C  Fba  U  {q}.  The  chain  D  must  contain 

H (Z).  Let  h'z  be  the  image  of  hx  under  the  mapping  point  q  or  else  C  would  be  the  boundary  of  a  chain  in  Fba 

H(X')  -3  H(Z').  Element  hx  is  the  image  of  both  hz  and  h  would  be  0.  Since  all  points  in  D  other  than  q  lie 

and  h'z  under  the  respective  mappings  H(Z)  -3  II (X )  in  F^’  ,  point  q  is  in  cl(Cj').  Since  D  connects  C  and  q  in 

and  H(Z')  -3  II (X).  cl (F_’a ) ,  cycle  C  and  point  q  lie  in  the  same  component 

By  Mayer- Vietoris,  the  sequence  ofcl(F*).  □ 


H(Y)  -3  H(Z)  ffi  H(Z')  -3  H(X) 

Theorem  1.  Let  f,g  :  fl  — »  R  be  continuous,  point 
is  exact.  Since  the  mapping  H(Z)  ®  H(Z')  -3  H(  X)  destructible,  LR  functions  on  fl  where  \  f  —  g\  <  8.  If  q 

sends  ( hz  ffi  h'z)  to  ( hx  —  hx)  =0,  element  ( hz  ffi  h'z)  destroys  non-zero  hf  G  i?(/_1(a))  and  \f(q)  —  a\  >  28 
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Figure  4:  Sets  Fba.  X  =  {x  £  tt  :  a  +  S  <  g(x)  and 
fix)  <b},  OH and  GH;+S. 

and  a  is  the  connected  component  of  c1(-F{^+2<5)  con¬ 
taining  q,  then  a  contains  a  point  q'  for  g  which  destroys 
some  hg  £  H{g^1{a  +  5))  and  f(q)  —  5  <  g{q')  < 
f(q)+S. 

Proof.  Let  b  equal  f{q).  Let 

I  =  {ie!!:a  +  i5<  g(x)andf(x)  <  b}. 

Since  \f(x)  —  g{x) \  <  S  and  b  —  a  >  25,  space  X  is 
properly  bounded.  Note  that  j  C  A'  C  F_ba  and 
so  there  exist  homomorphisms  H  (G^+Ss  )  -»•  H(X)  -3 
H(F_ba)  induced  by  inclusions. 

The  following  commutative  diagram  gives  the  rele¬ 
vant  mappings  between  homology  groups: 

hg  £  hf  £ 

Hig-Ha  +  S))  H(f-\a )) 


H(Gba+s5) - Jh\x) - 3  h'f 


K  £  H{G%+5)  H(X  U  {<?})  —  H(Fba  U  M) 


H(G™d) 

The  value  b'  will  be  defined  below. 

Let  h'f  £  H  ( F_a )  be  the  image  of  hf  under  the  map¬ 
ping  H(f-L(a))  -3  H(Fba).  Since  h'f  is  destroyed 
at  /_1(6),  element  h'f  is  non-zero.  By  Lemma  4,  el¬ 
ement  h'f  is  the  image  of  some  element  hx  £  if  (A) 
which  is  destroyed  by  q  and  is  the  image  of  some  hg  £ 


H(g^1(a  +  5)).  Since  h'f  is  non-zero,  elements  hx  and 
hg  are  non-zero. 

Since  hx  is  destroyed  by  q ,  the  mapping  if  (A)  -3 
H{ X  U  {</})  -3  H (Gff+s)  sends  hx  to  zero.  Thus  the 
composition  of  mappings  H(g^1(a  +  S))  -3  H(X)  -3 
H (G'tf+d)  sends  hg  to  hx  to  zero. 

By  Corollary  3,  there  exists  a  point  q1  £  Q  such  that 
hg  is  destroyed  by  q1  (i.e.,  the  image  of  hg  under  the 

mapping  H(g^1(a  +  6))  -3  H{Ga+s)  is  destroyed  by 
q' .)  Let  b'  equal  f(q').  Since  is  a  subset  of  X, 

the  image  of  hg  under  the  mapping  H{g  l{a  +  5))  -3 
H(Gba+Ss)  is  non-zero  and  so  b'  >  b  —  d.  Let  h'g  be 
the  image  of  hg  under  the  mapping  H{g  l{a  3-  5))  -3 

HiGi+s). 

Element  hg  is  generated  by  some  cycle  C  in  g  1  (a  + 
S).  Since  the  image  of  hg  is  hx  £  II (X)  under  the  map¬ 
ping  H(g  ^(a  +  5)  — ¥  H(X),  cycle  C  also  generates 
hx.  Similarly,  cycle  C  generates  h'g  £  H(Gba+s)  and 
h'f£H(Fba). 

Since  hx  is  destroyed  by  q ,  cycle  C  is  the  boundary 
of  some  chain  D  C  X  U  /_1(6).  Since  \f  —  g\  <  6, 
set  X  U  /_1(6)  is  a  subset  of  G^g  U  g^1^  +  S)  and 
so  chain  D  is  a  subset  of  G^j^  U  g^x{b  +  5).  Since 
h'g  €  H(G^+S)  is  non-zero,  space  G^-  U  g~l(b  +  5) 
cannot  be  a  supspace  of  G*+5 .  Thus  b'  is  at  most  b  +  5. 

As  noted  above,  cycle  C  generates  h'g  £  H(G^+S) 
and  h'f  £  H{F^f).  Point  q  destroys  h'f  which  is  gener¬ 
ated  by  C.  By  Lemma  6,  point  q  must  lie  in  the  con¬ 
nected  component  of  cl(F_a)  containing  C.  Similarly, 
since  q'  destroys  h'g,  point  q'  must  lie  in  the  connected 
component  of  cl(G*+^)  containing  G.  Since  cl(F*)  C 
cl (Fba+2d)  and  cl(G60'+5)  C  cl(G^)  C  c\(Fb+2S), 
points  q  and  q'  must  lie  in  the  same  connected  compo¬ 
nent  o  of  cl(F*+2<5).  □ 

5  Computing  persistence 

Theorem  1  can  be  used  to  compare  two  real  valued  func¬ 
tions  /  and  g  defined  on  a  topological  space  Q.  The  key 
computation  to  apply  Theorem  1  is: 

(i)  determine  if  a  point  p  which  destroys  some  h  £ 
H (f-1  (a))  has  persistence  greater  than  7. 

We  use  Betti  numbers  and  their  persistent  counter¬ 
parts  to  compute  (i).  Recall  that,  for  a  homology  group 
H (A'),  the  Betti  number  is  /3(A)  =  dim  H( X).  It  gives 
the  number  of  generators  in  H( A).  The  persistent  Betti 
numbers  relate  the  homology  classes  of  one  space  into 
the  other.  For  Y  C  A,  let  H ^  be  the  image  of  the  map 
H(Y )  -3  if  (A)  induced  by  inclusion  Y  -3  X.  De¬ 
fine  P(Y,X)  =  dim  if  y.  In  words,  8 (Y.  X )  counts  the 
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number  of  non-zero  generators  of  H  (Y)  that  remain  so 
in  the  larger  space  X. 

We  discuss  the  computations  for  the  function  /.  It  is 
clear  that  similar  computations  are  needed  for  g  as  well. 
In  general,  for  a  point  p  and  a  value  a  <  b  =  f(p) 
we  want  to  compute  if  an  element  of  H(f~1(a))  gets 
destroyed  by  p.  Let 

=  /?(/  1  (a), Kba)  and 

K  =  ^r1W,^p)uM). 

Note  that  p  is  a  point  whereas  a  and  b  are  real  val¬ 
ues.  The  number  Bb  counts  the  number  of  generators 
of  H(f^1(a))  surviving  in  F b  and  counts  the  num¬ 
ber  of  generators  of  H(f~1(a))  surviving  in  Fb  U  {p}. 
Therefore, 

<  =  &-  K 

counts  the  number  of  generators  of  H(f~1(a))  de¬ 
stroyed  by  p.  So,  if  tt?  >  0,  we  have  a  generator  of 
H (/_1  (a))  that  is  destroyed  by  p  where  /(p)  =  b. 

Notice  that  Theorem  1  can  be  applied  to  any  point 
p  £  0  and  any  value  a  <  / (p) .  However,  for  a  canonical 
computation  one  can  focus  on  the  critical  points  of  the 
functions.  We  define  the  critical  points  of  /  as  follows. 

Definition  1.  A  point  p  £  0  is  critical  for  f :  f  1  — >  K  if 
H  (Fb)  — »  H  (Fb  U  {p})  is  not  an  isomorphism  for  some 
a,  b  £  R  where  either  f(p)  =  a  or  f(p )  =  b. 

The  above  definition  of  critical  points  is  similar  to  that 
of  Cohen-Steiner  et  al.  [3]  with  a  distinction  that  the  rel¬ 
evant  space  is  an  interval  set  whose  lower  level  may  not 
be  at  —oo.  Also,  notice  that  the  space  Fb  could  be  above 
or  below  the  level  of  /(p). 

Compute  the  critical  points  of  /  according  to  the 
above  definition.  Let  po,Pt--Pk  be  the  critical  points 
of  /  ordered  according  to  the  increasing  values,  that 
is,  f(pi)  >  f(pi-i)  for  all  i  >  0.  We  compute  the 
persistence  for  these  critical  points  p*  as  fol¬ 

lows.  For  1  <  i  <  A:  —  1,  let  Oj  be  a  value  with 
f(p%- i)  <  at  <  f(pi).  Compute  7Ta-  for  any  pair  i,j 
where  j  >  i  >  0.  Since  Ttvf  is  constant  for  all  a  where 
/ (pi  i )  <  a  <  f(pi),  if  7tai  is  greater  than  0,  then  the 
persistence  Tlf(pj)  is  at  least  f{pj)  —  f(j>i-i)-  Thus  we 
compute  II f(j>j)  as 

max  ~  f(Pi- 1)|  so  that  tt«  >  0. 

Similarly,  we  can  compute  the  critical  points 
qo,  qi, qn  and  a  set  of  intermittent  values 
bi,b^n..,bn-i  for  the  function  g.  The  persistence 
of  a  critical  point  q  of  g  is  measured  similarly  by  IIa(g). 

To  compare  /  and  g ,  we  check  if  any  critical  point 
p  of  f  has  II /  (p)  greater  than  a  user  supplied  parame¬ 
ter  r.  If  so,  we  search  for  a  critical  point  q  of  g  in  the 


connected  component  of  cl(F^^j+^)  so  that  IIs(q)  >  r 
and  |  /(p)  —  g(q)\  <  j.  If  r  >  28 ,  such  a  q  exists  by 
Theorem  1 . 

5.1  PL  case 

Assume  there  is  some  finite  triangulation  of  !!  such  that 
/  and  g  are  linear  on  each  simplex  of  the  triangulation. 
Functions  /  and  g  are  LR  (locally  retractible),  but  not 
necessarily  point  destructible.  The  critical  points  of  / 
and  g  are  located  at  the  triangulation  vertices.  A  small 
perturbation  of  the  scalar  value  at  each  triangulation  ver¬ 
tex  and  the  linear  interpolation  of  those  values  over  the 
triangulation  simplices,  gives  new  piecewise  linear  func¬ 
tions  which  are  point  destructible. 

Carlsson  and  Zomorodian  in  [9]  show  how  to  compute 
persistent  Betti  numbers  for  homology  groups  of  filtered 
simplicial  complexes  over  any  field.  However,  spaces 
F_ba  and  F_ba  U  {p}  are  not  closed.  To  compute  their  per¬ 
sistent  Betti  numbers,  Bb  and  \va,  we  collapse  them  onto 
closed  sets  which  are  simplicial  complexes. 

Consider  the  space  F_ba .  Let  t  C  Q  be  a  simplex  with 
a  point  in  this  space.  If  each  vertex  v  of  t  either  has 
f(v)  >  b  or  a  <  f(v)  <  b,  the  subset  t  FI  Fba  can 
be  collapsed  to  the  face  of  t  made  by  the  vertices  whose 
values  lie  in  [a,  b).  This  cannot  be  done  for  simplices  that 
cut  across  the  levels  of  a  and  b.  These  simplices  have 
vertices  with  values  above  b  and  also  below  a.  For  such 
a  simplex  we  take  an  edge  e  =  {u,v}  where  f(u )  < 
a  and  f(v )  >  b  and  consider  a  point  x  on  this  edge 
where  a  <  f(x)  <  b.  We  divide  t  by  starring  from  x 
to  all  its  vertices.  After  subdividing  all  such  simplices 
we  obtain  a  subdivision  H  of  fl  which  has  no  simplex 
cutting  across  the  interval  [a,  b].  Consider  the  simplicial 
complex  made  by  the  collection  of  simplices  in  O  that 
have  all  vertices  with  values  in  [a,b).  The  underlying 
space  of  this  simplicial  complex  is  a  deformation  retract 
of  F_ba  and  therefore  has  homology  groups  isomorphic  to 
that  of  F_a  ■ 

6  Maxima 

We  show  that  the  neighborhoods  of  local  maxima  with 
large  persistence  are  pairwise  disjoint.  This  enables  us 
to  establish  a  matching  of  such  critical  points. 

The  idea  of  the  proof  is  as  follows.  Consider  two  local 
maxima,  p,p'  £  fl,  where  f(p)  <  f(p')-  If  p  destroys 
non-zero  h  £  Hk  (/_1  (a)),  then  k  equals  d  —  1.  Let  a  be 
the  connected  component  of  Fff  containing  p.  If  a  is  a 
manifold  with  boundary,  then  U  {p}  must  contain 
a.  Since  f(jp)  <  f(p'),  set  U  {p}  does  not  contain 
p'  and  therefore  point  p'  is  not  in  a. 
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Lemma  7.  Let  Fff  be  an  oriented  d-manifold  with 
boundary.  If  p  is  a  local  maximum  and  p  destroys  non¬ 
zero  h  £  Hk  ),  then  k  equals  d  —  1. 

Proof.  See  appendix.  □ 

As  previously  noted,  we  always  mean  path  connected 
when  referring  to  connected  components. 

Lemma  8.  Let  M  be  a  connected,  oriented  d-manifold 
with  non-empty  boundary.  Let  D±,D 2, . . .  ,Dk  be  the 
connected  components  ofdM  with  orientation  inherited 
from  M.  If  a\D\  +  a2D2  +  . . .  akDk  generate  the  zero 
element  of  Ha-i(M),  then  01  =02  =  . . .  =  a*. 

Proof.  The  sequence  Hd(M)  — >  Hd(M,dM)  — » 
Ha-i(dM)  — >  Ha-i(M)  is  exact  [6,  Theorem  2.16, 
p.  117],  Since  the  boundary  of  M  is  not  empty,  the 
homology  group  Hd(M)  is  zero.  The  homology  group 
of  Hd(M,dM)  is  G,  the  ground  ring  of  the  homology 
group.  The  map  Hd(M,dM)  — >  Hd~i(dM)  is  the  con¬ 
necting  homomorphism.  It  maps  G  to  h  £  Hd-i(dM) 
which  is  generated  by  D1  +  D2  +  . . .  +  Dk.  Since  the 
mapping  is  exact,  the  image  of  G h  under  the  mapping 
Hd-i(dM)  is  zero.  Moreover,  only  el¬ 

ements  in  Gh  map  to  zero.  Thus,  if  aiDi  +  a2D2  + 
. . .  akDk  generate  the  zero  element  of  Hd-i(M),  then 
a  1 1)  1  +  a2D2  +  . . .  a/, I) i-  must  generate  an  element  of 
G h  in  Hd-i(dM)  and  so  oq  =  a2  =  . . .  =  ak-  □ 

Lemma  9.  Let  M  be  a  connected,  oriented  d-manifold 
with  non-empty  boundary.  If  dM  C  M1  C  M 
and  Hd-i(dM)  — *■  Hd-i(M')  takes  non-zero  h  £ 
Hd~i(dM)  to  zero,  then  M'  equals  M. 

Proof.  Assume  M'  does  not  equal  M.  Let  p  he  a  point 
in  M  —  M'.  Let  B  be  an  open  topological  ball  containing 
p  whose  closure  does  not  intersect  dM.  There  exists  a 
deformation  retract  from  M  —  {p}  to  M  —  B.  Thus 
H(M  —  {p})  is  isomorphic  to  H(M  —  B). 

The  mapping  Hd-i(dM)  Hd  i(M')  -y 

Hd-i(M  —  {p})  — Hd~i(M  —  B)  sends  h  to  zero 
in  —  B).  Let  ho  be  the  element  of  the  homol¬ 

ogy  group  of  Hd-i(dM)  generated  by  dM  with  orienta¬ 
tion  inherited  from  M.  By  Lemma  8,  element  h  equals 
a  ho  for  some  non-zero  a.  Let  h  u  be  the  element  of 
Hd-\{M  —  B)  generated  by  OB  with  orientation  inher¬ 
ited  from  M  —  B.  Let  h'  be  the  image  of  h  under  the 
map  Hd-i{dM)  — >  Hd-i{dM  U  dB).  The  element 
h  and  hence  h'  is  sent  to  zero  in  Hd-i(M  —  B).  By 
Lemma  8,  element  h'  equals  /3(ho  +  hs)  for  some  non¬ 
zero  B.  Thus  aho  equals  fi(ho  +  he)-  Since  ho  and  hs 
are  linearly  independent,  a  and  B  are  both  zero  implying 
h  is  a  zero  element,  a  contradiction.  It  follows  that  M 
equals  M' .  □ 


Let  a?  (7)  represent  the  connected  component  of 
—Tip)- 7  containing  P-  We  prove  that  the  neighborhoods 
of  (7)  of  points  with  persistence  greater  than  7  are  pair¬ 
wise  disjoint.  (See  Figure  5.) 

Theorem  2.  Let  f  :  0  — »  R  be  a  continuous  function 
such  that  Fff  is  a  (d  —  l)-manifold  with  boundary  for 
all  but  a  finite  number  of  a.  If  points  po,Pi  €E  fl  are  lo¬ 
cal  maxima  with  persistence  greater  than  7,  then  ojJ0  (7) 
does  not  intersect  o^  (7). 

Proof.  Let  po-  Pi  G  fi  be  local  maxima  with  persistence 
70, 71,  both  greater  than  7.  Without  loss  of  generality, 

assume  that /(p0)  <  fipi). 

Assume  that  ofQ  (7)  intersects  o^  (7).  Since  is  a 
(d  —  1) -manifold  for  all  but  a  finite  number  of  a ,  there 
is  some  7'  >  7  such  that  j0  >  l'  and  ^  >  7'  and 
—T(po)  7'  is  a  (d  ~  l)-manifold  with  boundary.  Since 
of 0  (7)  intersects  o £  (7),  set  ofQ  (7')  intersects  o^  (7'). 

Since  /(p0)  <  /(pi),  set  Ff{po)_y  contains 

Ejfa)--/-  ■  Since  (7')  intersects  (7'),  set  a/0  (7') 
contains  ofpi  (7').  Thus  ofQ  (7')  contains  pi. 

By  Lemma  4,  point  p0  destroys  some  non-zero  ele¬ 
ment  of  Hk(f~1(f(po)  —  7')).  By  Lemma  7,  k  equals 
d—  1.  By  Lemma  10  (appendix),  point po  destroys  anon- 
zero  elements  of  Hd~i(dof0( 7')).  In  Lemma  9  putting 
M  =  ofQ  (Y)  and  M'  equal  the  connected  component  of 

— Z(po)-7'  containing p0,  we  conclude  F_fj%°0 j_7<  U  {p0} 
contains  ofQ  (7')  and  thus  contains  p\.  However,  since 

fipo)  <  /(Pi)-  set  E£I'po'>  U  {po}  does  not  contain  pi. 
Thus,  ofQ  (7)  does  not  intersect  o^  (7).  □ 

Our  final  theorem  gives  relationships  between  neigh¬ 
borhoods  of  local  maxima  of  /  and  of  g. 

Theorem  3.  Let  /,  g  :  fl  — >  R  be  continuous  functions 
such  that  Fff  and  Gff  are  (d—1) -manifolds  with  bound¬ 
ary  for  all  but  a  finite  number  of  a  and  \f  —  g\  <6.  Let 
p  £  17  be  a  local  maxima  of  f  and  let  q  and  q1  be  local 
maxima  of  g  such  that  p,  q ,  q'  have  persistence  greater 
than  7  and  |/(p)  —  g(q) |  <  6  and  |/(p)  —  g(q')\  < 

(i)  If  o.f  (7  —  26)  intersects  cr|( 7  —  25),  then  <7®  (7) 
contains  of  (7  —  2$). 

( ii )  If  of  (7  —  26)  intersects  ct|  (7  —  25),  then  of  (7  —  25) 
does  not  intersect  o9q,  (7  —  25). 

Proof  of  (i).  Let  y  be  a  point  in  of  (7  —  25)  0  a®  (7  —  2<5) 
and  2  be  any  point  in  of  (7  —  25).  Set  of  (7  —  25)  is  path 
connected,  so  there  is  a  path  (  C  of  (7  —  25)  from  y  to  2. 
Since  (  C  <r|(7  -  25),  f(x)  >  f(p)  —  7  +  25  for  every 
point  x  £  C-  Since  |  f{x)  —  g(x)  \  <  5  for  all  x  £  U  and 
\f(p)  —  g(q)\  <  5,  it  follows  that  g  (x)  >  g  ( q )  —  7  for  all 


x  £  Q.  Thus  Q  C  a9  (7)  and  z  lies  in  a^('y).  This  holds 
for  all  z  G  a 7  -  25)  so  <r|( 7)  contains  <7^(7  -  25). 
(See  the  neighborhoods  of  P2  and  q-2  in  Figure  5.)  □ 

Proof  of  (ii).  By  Theorem  2,  o9(  7)  and  a9,  (7)  are  dis¬ 
joint.  By  (i)  above,  if  <7^  (7  —  25)  intersects  o9  (7  —  25), 
then  cr|  (7)  contains  <7^  (7  —  25) .  Similarly,  if  (7  —  2 5) 
intersects  <7®,  (7  —  25),  then  <7®,  (7)  contains  cr?  (7  —  25). 
However,  ct|  (7)  and  o9 ,  (7)  are  disjoint.  Thus,  if  (7  — 
25)  intersects  (7  —  25),  then  <7^  (7  —  25)  does  not  in¬ 
tersect  <7®  (7  —  25).  (See  Figure  5.)  □ 

7  Matching 

We  assume  that  /,  g  :  fi  — >  R  are  continuous,  point 
destructible,  LR  functions  such  that  Fff  and  Gff  are 
(d  —  1  )-manifolds  with  boundary  for  all  but  a  finite  num¬ 
ber  of  a.  Let  Mf  and  Mg  be  the  set  of  local  max¬ 
ima  of  /  and  g,  respectively,  and  let  Mf( 7)  C  Mf 
and  Mg( 7)  C  Mg  be  the  set  of  local  maxima  of  / 
and  g,  respectively,  with  persistence  greater  than  7.  We 
would  like  to  match  points  in  Mf  (7)  with  close  points  in 
Mg  (7)  in  the  sense  of  Theorem  1 .  However,  there  may 
be  no  such  matching.  In  fact,  /  may  contain  a  set  of 
critical  points  with  persistence  a  little  bit  above  7  while 
nearby  critical  points  in  g  all  have  persistence  a  bit  be¬ 
low  7.  Thus,  Mf( 7)  can  contain  any  number  of  points 
while  Mg(  7)  is  empty!  Instead  of  matching  Mf(  7)  and 
Mg( 7)  only  with  each  other,  we  allow  them  to  match 
with  points  with  slightly  less  persistence. 

We  say  that  a  partial  matching  of  M f  with  Mg  covers 
Mf( 7)  and  Mg( 7)  if  all  points  in  Mf  (7)  and  Mg( 7)  are 
matched.  As  before,  let  <7^(7)  and  a9  (7)  be  the  con¬ 
nected  components  of  Eqpp)  7  ar,d  —9(9)- 7  containing 
p  and  5,  respectively.  A  partial  matching  of  Mf  and  M3 
is  (ct,  j3)-close  if  for  each  pair  ( p ,  </)  where  p  £  Mf  and 

£  Mg,  point  q  lies  in  a? (a)  and  point  p  lies  in  <7 9  (a) 
and  |  f  {p)  —  g(q)\  <  3. 

Assume  that  \f(x)  —  g(x) \  <  7/4  for  all  x  €  fi.  We 
will  find  a  partial  matching  of  Mf  with  Mg  which  covers 
Mf( 7)  and  Mg( 7)  and  is  (7, 7/4)-close. 

The  algorithm  is  as  follows.  For  each  point  p  £ 
Mf  (7),  we  compute  ap  =  ^(7/2).  Similarly,  for 
each  £  Mg (7),  we  compute  <7g  =  a|( 7/2).  By 
Theorem  3,  each  <7p  intersects  at  most  one  aq  where 
|/(p)  ~ 9(l)\  <7/4  and  vice  versa.  If  op  intersects  such 
a  oq ,  then  match  p  with  q.  If  not,  then  match  p  with  some 
q '  G  Mg( 7/2)  lying  in  such  that  |/(p)  -  /(g')l  < 
7/4.  (By  Theorem  1  such  a  q1  exists.)  Similarly,  if  aq 
does  not  intersect  any  op,  match  q  with  p1  £  Mf( 7/2) 
lying  in  aq  such  that  \f(q’)  -  f{p) \  <  7/4. 

We  claim  that  algorithm  MatchPersistentMax 
matches  all  maxima  with  persistence  more  than  7: 


MatchPersistentMax!  H,  /,  7) 

/*/,£/:  fi  ->  M  */ 

1  Compute  sets  Mf  (7)  and  Ms  (7); 

2  For  each  point  p  £  Mf( 7),  compute  <7P  = 

(7/2); 

3  For  each  point  g  £  Ms( 7),  compute  = 
°g(  7/2); 

4  For  each  point  p  £  Mf( 7)  and  5  G  Ma( 7),  if  FI 
<7g  7^  0  and  |  f(p)~  d(q)\  <  7/4,  then  match  p  with 

5  For  each  unmatched  p  £  Mf  (7),  match  p  with  </'  G 
M9(7/2)  n  where  |/(p)  -  /(g')|  <  7/4; 

6  For  each  unmatched  5  G  Ms( 7),  match  5  with  p'  G 

Mf{ 7/2)  n  where  |5(9)  -  g(p')|  <  7/4- 


Figure  5:  Local  maxima  and  their  neighborhoods. 
Solid  lines  around  points  p;-  and  q j  are  neighborhoods 
aPi  (7/2)  and  t7!;  (7/2)-  Dotted  lines  are  neighborhoods 
api('y)  and  a|i(7)-  Neighborhood  aj/2(7/2)  intersects 
af2  (7/2)  so  P2  matches  with  52-  Point  pi  matches  with 
some  point  (not  shown)  from  Ma( 7/2)  and  points  gi  and 
^3  match  with  points  (not  shown)  from  Mf( 7/2). 

Proposition  1.  If  \f  —  g\  <  7/4,  then 

MatchPersistentMax(Q,  /,p,7)  produces  a 
partial  matching  of  Mf  with  Mg  which  covers  M f  (7) 
and  Mg  (7)  m/c/z  f/znf  every  matched  pair  (p,q)  where 
p  £  Mf  and  q  £  Mg  is  [pf,'y/A)-close. 

Proof  By  Theorem  2,  the  ap  =  <7^ (7/2),  p  G  Mf( 7), 
are  pairwise  disjoint  and  the  <t,j  =  a|( 7/2),  5  G  Mg( 7) 
are  pairwise  disjoint.  By  Theorem  3,  each  ctp  intersects 
at  most  one  oq  and  vice  versa.  Thus  Step  4  gives  a  one 
to  one  partial  matching. 

By  Theorem  1,  op  contains  some  point  q1  £  M( 7/2) 
such  that  |/(p)  —  g{q')\  <  7/4.  Since  (p,q')  is  not 
matched  in  Step  4,  point  (/'  is  not  in  Mg( 7).  Thus  point 
<f  is  not  matched  in  Step  4.  Since  does  not  intersect 
any  ap' ,  p'  G  Mf( 7),  point  g'  is  matched  to  at  most  one 
p  in  Step  5.  Similarly,  point  p'  in  Step  6  is  not  matched  in 
Steps  4  and  5  and  is  matched  to  at  most  one  q  in  Step  6. 
Thus  the  matching  is  one  to  one  and  covers  all  of  My  (7) 


9 


and  Mg  (7). 

It  remains  to  show  that  for  each  match  ( p ,  q),  set 
ap  (l)  contains  q  and  cr|  (7)  contains  p.  If  p  and  q  are 
matched  in  Step  4,  then  a jj  (7/2)  intersects  and  a9  (7/2). 
This  holds  true  even  if  p  and  q  are  matched  in  Steps  5 
or  6.  By  Theorem  3  with  S  =  7/4,  <7^(7)  contains 
<7®  (7/2)  which  contains  q  and  a9( 7)  contains  ajj( 7/2) 
which  contains  p.  Since  points  p  €  Mf  and  q  €  Mg  are 
only  matched  if  \f(p)  —  g(o)\  <  7/4,  the  matching  is 
(7,7/4)-close.  □ 

8  Discussions 

Results  on  stability  of  topological  persistence  can  be 
used  in  shape  matching.  If  we  take  a  dense  point  sample 
from  the  boundary  of  a  shape,  the  distance  functions  to 
the  shape  boundary  and  and  its  point  sample  are  simi¬ 
lar.  Therefore,  if  we  have  two  similar  shapes,  the  dis¬ 
tance  functions  defined  by  their  point  samples  are  simi¬ 
lar.  As  observed  in  previous  works  [1,2],  the  results  on 
persistence  apply  to  such  functions.  Our  results  in  this 
paper  have  some  notable  connections  to  a  shape  match¬ 
ing  algorithm  proposed  by  Dey  et  al.  [4].  According  to 
our  results,  we  can  expect  that  similar  shapes  have  simi¬ 
lar  structures  for  maxima  with  large  persistence  in  terms 
of  the  interval  sets.  The  algorithm  in  [4]  uses  maxima 
and  their  stable  manifolds  for  matching.  We  suspect  that 
these  stable  manifolds  are  playing  the  role  of  connected 
components  as  suggested  in  this  paper.  Perhaps  the  per¬ 
formance  of  the  matching  algorithm  in  [4]  now  can  be 
improved  and  better  explained  by  our  results.  We  plan 
to  address  this  issue  in  future  work. 

Persistence  diagrams  [3]  can  be  used  to  match  shapes 
based  on  the  critical  values  of  the  distance  function. 
Does  adding  critical  points  increase  the  discrimination 
of  this  matching?  Figure  6  provides  such  an  example. 
The  two  unsimilar  shapes  in  this  figure  has  matching  per¬ 
sistence  diagrams.  The  distance  function  for  each  shape 
contains  three  local  maxima  and  two  saddle  points.  For 
each  shape,  two  of  the  local  maxima  and  one  of  the 
saddle  points  has  long  persistence  while  the  other  local 
maxima  and  saddle  point  have  short  persistence.  (Persis¬ 
tence  of  saddle  points  which  “create  homology”  groups 
is  defined  in  [3].)  Thus  the  persistence  diagrams  match. 

Now  assume  that  the  two  shapes  are  registered  so  that 
each  qi  lies  on  top  of  77.  Even  with  such  a  registration, 
the  local  maxima  77  does  not  match  with  <7  since  77 
and  (p  have  different  persistence  and  77  does  not  match 
with  q-2  or  q:>  since  they  are  not  in  the  neighborhood 
of  pi  in  the  domain.  Thus,  our  matching  based  on 
closeness  of  local  maxima  in  both  the  range  and  domain 
distinguishes  these  shapes. 


Figure  6:  Two  shapes  with  matching  persistence  dia¬ 
gram.  Matching  based  on  closeness  in  both  the  range 
and  domain  distinguishes  these  shapes.  Local  maxima 
77,  P3 •  Qi,  93  and  saddle  points  p5,  q5  have  long  per¬ 
sistence.  Local  maxima  p2,  qi  and  saddle  points  7)4,  q.i 
have  short  persistence. 
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9  Appendix 

Proof  of  Lemma  1.  Assume  function  /  :  PI  ->  H  is 
LR.  Choose  ei  >  0  such  that  f^1(a)  is  a  strong  de¬ 
formation  retract  of  {x  :  a  —  e  <  f(x )  <  a  +  e} 
for  all  e  <  ei  and  a  +  e  <  b.  Similarly,  choose 
C2  >  0  such  that  /  l(b)  is  a  strong  deformation  retract 
of  {x  :  b— e  <  /(x)  <  fr+e}  for  all  e  <  e2  and  a  <  b—c. 
Let  €q  be  the  minimum  of  ei  and  e->  and  b  —  a. 

For  any  e  <  eo-  let  <Pc  t>e  the  mapping  from  {a  :  a  — 
e  <  f(x)  <  a  +  e}  x  I  to  f^1(a)  representing  the 
strong  deformation  retract  of  {a  :  a  —  e  <  x  <  a  +  e}  to 
f^1(a).  Let  <j>\  be  the  mapping  from  {6  :  b— e  <  f(x)  < 
b  +  e\  x  I  to  /  1  (b)  representing  the  strong  deformation 
retract  of  {6  :  b  —  e  <  x  <  b  +  e}  to  /_1(6).  Define 

{</>“(x,i)  for  a  —  e  <  x  <  o, 
x  for  a  <  x  <  b, 

<j>b(x,t)  for  b<x<b  +  e. 

<f>t  is  constant  on  {x  :  a  <  x  <  b}  and  continously 
deforms  {x  :  a  —  e  <  x  <  a}  and  {x  :  b  <  x  <  b  +  e} 
onto  {x  :  a  <  x  <  b}.  Thus  {x  :  a  <  f(x)  <  6}  is  a 
strong  deformation  retract  of  {x  :  a  <  f(x)  <  5  +  e} 
and  {x  :  a  <  f(x)  <  6}  is  a  strong  deformation  retract 
of  {x  :  a  —  e  <  f(x)  <  b}.  □ 

Proof  of  Lemma  2.  Since  H(f^1(a))  -»  H{F_ “) 

sends  h  to  zero,  element  h  is  the  boundary  of  some  chain 
C  C  Fff.  Chain  C  is  compact.  (See  [8,  p.  71].)  Thus 
{f(x)  :  x  G  C)  is  compact  and  has  a  maximum  value 
b'.  Since  C  C  F% ,  the  mapping  H(f-1(a))  ->  H(F*) 
sends  h  to  zero. 

Let  b  equal  liminf{&  :  H(f~1(a))  — >  H(F^a)  sends 
h  to  zero}.  Note  that  b  <  b'.  Since  /  is  LR,  H(f~l(a)) 
is  isomorphic  to  H(Ffffl)  for  sufficiently  small  e  and 
thus  H(f~1(a))  — »  iT(F“+e)  does  not  send  h  to  zero. 
Thus  b  is  strictly  greater  than  a. 

Let  h'  be  the  image  of  h  under  the  mapping 
H (/_1  (a))  — »  H (F_ba).  If  h'  were  zero,  then  h  would  be 
the  boundary  of  some  chain  C'  C  F_ba.  Since  C'  is  com¬ 
pact,  chain  C'  would  also  be  a  subset  of  F_ba  for  some 
b  <  b,  contradicting  the  choice  of  b.  Thus  h'  is  non¬ 
zero. 

We  show  that  h'  is  destroyed  by  /_1(6).  Since 
/  is  LR,  there  is  some  eo  >  0  such  that  H(F_ba  U 
f^b(b))  is  isomorphic  to  H(Fba]~e )  for  all  e  <  e0- 
If  H(Fba)  ©>  H(Fba  U  f-^b))  does  not  map  h'  to 
zero,  then  H(Fba )  — >  H(F_ba+e)  does  not  map  h'  to 
zero  for  all  e  <  eo,  and  b  does  not  equal  liminf{6  : 
H(f~l(b))  — >  H(F^)  sends  h  to  zero}.  Thus, 

H(F ]£)  — >  H(F_ba  U  /“*(&))  maps  h'  to  zero  and  h'  is 
destroyed  by /_1  (6).  □ 


Proof  of  Lemma  7.  Since  Fff  is  an  oriented  d- 
manifold  and  p  is  a  local  maximum,  some  neighbor¬ 
hood  Np  of  p  is  homeomorphic  to  Rrf  and  all  points  in 
Np  —  {p}  have  value  less  than  f(p).  Let  B  be  the  unit 
ball  in  3,\  and  let  Bp  be  its  image  under  the  homeomor- 
phisrn  from  to  Np.  Since  all  points  in  Np  —  {p}  have 
value  less  than  /(p),  they  are  all  in  F_^p\ 

By  the  Mayer- Vietoris  Theorem,  the  sequence 

Hk(Bp-{p})  -j.  Hk(Bp)®Hk(Ei^)  -j.  Hk(F^U{p}) 

is  exact.  Since  the  mapping  — >  Hk(F ^  U 

{p})  sends  h  to  zero,  the  mapping  Hk(Bp)  © 
Hk{Fl (p))  -»■  Hk(F{^  U  {p})  sends  (0  ©  h)  to  zero. 
Since  the  sequence  is  exact,  element  (0  ©  h)  is  the  im¬ 
age  of  some  non-zero  h'  G  Hk(Bp )  under  the  map¬ 
ping  Hk(Bp  -  {p})  Hk(Bp )  ©  Hk{F{ W).  Since 
Hk(Bp  —  {p})  is  the  zero  group,  for  all  k  ^  d  —  1,  ele¬ 
ment  h1  must  be  in  II,  i  \(Bp  —  {p}).  Therefore,  h  is  an 
element  of  and  so  k  equals  d  —  1.  □ 

Lemma  10.  Let  DC  [V  he  topological  spaces,  let 
a[.  . . . ,  o'm  be  the  pathwise  connected  components  of 
Ll' ,  and  let  o  equal  for  each  i.  If  the  mapping 

Hk(Q)  Hk(Ll')  sends  non-zero  h  G  Hk(fl)  to  zero, 
then  there  exists  some  non-zero  hi  G  Hk(uf)  such  that 
the  mapping  Hk{af}  — >  Hk(Q')  sends  hi  to  zero.  More¬ 
over,  if  point  p  is  an  element  ofLl'  —  D,  and  the  mapping 
Hk(oi)  ->  Hk(Ll'  —  {p})  does  not  send  hi  to  zero,  then 
p  is  an  element  of  a[. 

Proof  Since  the  Oi  are  pairwise  disjoint,  the  homology 
group  Hk  (D)  is  isomorphic  to  Hk  (ay)  ©  ■  ■  ■  ©  Hk  (. am ). 
(See  [8,  Theorem  4.13,  p.  69].)  This  isomorphism  takes 
h  to  (hi  ©  ■  ■  ■  ©  hm)  where  hi  G  Hk(oi).  At  least  one 
of  these  hi  must  be  non-zero. 

Since  the  o(  are  pairwise  disjoint,  the  homology  group 
Hk(fl')  is  isomorphic  to  llk(a[  )©•••©  Hk(a'm).  The 
mapping  Hk  (0)  — »  Hk  (Pi1)  sends  h  to  zero,  so  the  map¬ 
ping 

Hk(oi)  ©  ■  ■  •  ffi  Hk(am)  — >  Hk(a[)  ©  ■  ■  ■  ©  Hk(a'm) 

sends  (hi  ©■ « *ffi/tm)  to  (0©-  ■  -fflO).  Thus,  the  mapping 
Hk(oi)  Hk(a[)  -©  iTfc(D'),  sends  non-zero  hi  G 
Hk(oi)  to  zero. 

Let  point  p  be  an  element  of  Pi' —PI  where  the  mapping 
Hk (<ji)  -»  Hk(Pl'  —  {p})  sends  hi  to  some  non-zero 
h!  G  Hk(Pl'  -  {p})-  If  p  k  a'i •  then  a[  C  PI'  —  {p}  and 
the  mapping  Hk(oi)  -©  Hk(a'i)  ©■  Hk(Pl'  -  {p})  sends 
hi  to  zero.  Thus  a\  contains  p.  □ 
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